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Abstract 

We study the associated production of J/ip +7 in fixed-target experiments with 
nuclear targets as well as at the relativistic heavy-ion collider (RHIC). We find that 
this process affords a very clean probe of p g , the ratio of gluon density in a heavy 
nucleus to that of a proton. The combined x range thus available can be used to 
discriminate between the predictions of different models of the EMC effect for p g . 
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The importance of a good knowledge of the parton densities in hadrons cannot be 
overemphasized in view of their significance for the physics that can be studied at both 
the pp and the relativistic heavy ion colliders. The parton distributions, which are often 
obtained by using nuclear targets, have a non-trivial dependence on the environment, 
viz. the well-known EMC effect [III |2J. While the nuclear effects on the quark densities 
can be directly studied in deep inelastic scattering (DIS) experiments, the gluon densi- 
ties g(x, Q 2 ) get determined only indirectly, through their effect on the QCD evolution 
of the structure function. A good determination of the gluon densities in nuclei is possi- 
ble only through studies of different hard scattering processes, dominated by gluons, in 
photon/lepton/hadron - nucleus collisions. A direct determination of the nuclear gluon 
densities will help in pinning down the correct model of the EMC effect. The theoretical 
basis of this effect still remains unclear as there are a large number of models 0, all 
of which are in reasonable conformity with the data available from DIS. However, these 
models give very distinct predictions for p g defined as, 

Pg = g A (x,Q 2 )/g^x,Q 2 ). (1) 

Thus experimental information on this ratio may help in clarifying the dynamical origin 
of the EMC effect. 

The standard processes that have been used to measure gluon densities in protons are 
large-pr jet production, heavy flavour production and direct photon/dimuon production 
at large in hadron-nucleon collisions, and lepto- and photoproduction of quarkonia with 
nucleon targets. The same processes can, of course, be used as an effective tool to determine 
the gluonic EMC ratio, p g [[3], ||. In fact the NMC collaboration Q has measured the ratio 
of gluon distributions in Sn and C, by studying inelastic J/ip production by scattering 200 
and 280 GeV muon beams off these targets. With the cuts imposed in this experiment, 
the values of x probed are in the range 0.05 to 0.15. It will be interesting to complement 
this information with gluon distributions from other processes mentioned above. With the 
exception of lepto- and photoproduction of J/ip, and jet production at lower p T values, 
all the other processes receive contributions from quarks as well. So one has to either 
choose special kinematic regions or find clever combinations of different cross-sections Q 
to separate the quark contributions from those of the gluon. 

In the context of the colour singlet model ||, it has been shown recently that 
a study of the associated production of J/ip and 7 in hadronic collisions affords a very 
clean determination of the gluon density. In this process, one looks for a photon recoiling 
against the J/ip produced at large pr- The authors of Ref. have studied this process for 
the pp and the HERA ep collider and their study has been later extended || to polarised 
pp scattering. In this note, we analyse the effectiveness of this process to probe the gluon 
densities in nuclei in fixed- target experiments and at the RHIC. We find that, owing to 
the lower values of y/s in fixed target experiments, the range of x values for which g A 
can be probed increases considerably. This can be done by studying the triple differential 
cross-section da/dpr dy n, dy J '^ for negative rapidities of the photon and J/ip. Moreover, 
when the J/ip is tagged through its leptonic decay the final state is particularly clean. 
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This further increases the efficacy of this process in probing g A at larger values of x, as 
compared with direct photon production ||. The study of this process at the RHIC will 
be shown to probe g A at the same small values of x as in lepto- and photoproduction of 
J/ip. Determination of p g using this process turns out to have two additional advantages: 
i) it is less sensitive to the theoretical uncertainties of the production mechanism, and ii) 
it is directly proportional to the ratio of experimentally measured cross-sections. 

Since some of the statements regarding the production process depend crucially on 
the hadronization mechanism in the production of J/ip, it is worthwhile summarising its 
main features. J / ip production is assumed to take place through the mechanism specified 
in the colour singlet model 0. In this model, one projects out the state with the correct 
spin, parity and colour assignments from the full cc production amplitude, to match the 
quantum numbers of the Jfij). This projection is done at the level of the hard scattering 
amplitude itself, which yields a multiplicative factor, R(0), the J/ip wavefunction at the 
origin in coordinate space. 

The basic subprocess, at the parton level, which is responsible for the final state with 
the J I ip and the photon at large px is 

g + g -> J/iP + 7. (2) 

The qq initial state does not contribute because in that case the produced cc pair is not 
a colour singlet. Thus this final state is accessible only to gluons, and therefore depends 
directly on the gluon distribution. Moreover, the P-state charmonia do not get produced 
in this process. This is so because the subprocess 

9 + 9 ^X^ + 1 (3) 
does not produce the J/ip and photon at large px, whereas the subprocess 

9 + 9 -» X + 7 (4) 

is disallowed because of C-invariance and colour neutrality of the final state. Thus, this 
process is purely gluon-initiated and one does not have to worry about x production and 
the subsequent decay of the x into a J/i>- It is worth noting that for the P-states it is 
the derivative of the bound-state wavefunction at the origin, and not the wavefunction 
itself, that is relevant. From simply this consideration, one expects x production to be 
suppressed by a factor of 50, given that the amplitudes for the processes in eqs. (0) and 
(^]) are of the same order of magnitude. Also, outside the framework of the colour singlet 
model the qq initial state can give rise to a final state containing a quarkonium and 7. 
However, this production mechanism will give rise to additional hadronic activity in the 
vicinity of the photon. Hence demanding a high-p^, isolated 7 will always ensure that the 
process is gluon-initiated. 

It may be noted that in the framework of the colour singlet model the formation of 
J/ip takes place on a perturbative time scale, being ~ 1/mj/^. This can have extremely 
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interesting implications for the J/ip suppression signal || of quark gluon plasma forma- 
tion, which may also have a comparable or even larger formation time. Thus the only 
effect the plasma can have in this case will be on the propagation of the J/ip. This 
makes the physics of confinement /deconfmement cease to have any significant effect on 
the experimentally observed suppression of the J/ip cross-sections. Instead the physics 
of its relative propagation through quark gluon plasma and a hadronic gas assumes an 
important role, bringing out further the model-dependence of the proposed signal. 

The colour singlet model detailed above has been used in the past to describe both 
leptoproduction and hadroproduction of J/ip, and is known to give a good description of 
the kinematical distributions || [L0| . However, it has been found that there is considerable 
uncertainty in the overall normalisation for J/ip production. This is probably due to the 
non-relativistic treatment of the J/ip in this model, as well as to the small mass of the 
charm quark. However, we propose to study ratios of the cross-sections for different targets 
and hence the effect of these K-factors will be reduced. Moreover, the data |J require a 
K-factor > 1.0. This would mean that our cross-sections can indeed be considered as a 
conservative estimate. 

The cross-section for the subprocess g + g — > J/ip + j is given in the colour singlet 
model as || |7j 



da _167iaa 2 s M,p\R(0)\' 
~dl ~ 27S 2 
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(5) 



where s, t and u are the Mandelstam variables for the subprocess, and the modulus 
squared, |-R(0)| 2 , of the wavefunction is related to the leptonic decay width by 

9M 2 

\R(0)\ 2 = T7r4 r (W - e + e~) = 0.544 GeV 3 , (6) 

if the value of the leptonic width of the J/ip is taken to be 5.36 keV. 

The triple differential cross-section da/dpr dipdy J l^ for the reaction B + A — >• J/ip + 
7 + X is then given by 

^^=2p T W(a; 1 )/(x 2 )|, (7) 

where g B (xi) and g A (x 2 ) are the gluon distributions in the beam B and the nuclear target 
A respectively, da/di is the subprocess cross-section given by eq. @ and y J ^ and y 1 are 
the rapidities, in the pA centre-of-mass frame, of the J/ip and the photon respectively. 
The momentum fractions x\ and 22 of the beam and the target carried by the gluons are 
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riven by the following kinematic relations: 



1 r_ v J/i> V 7i 

xi =^x T e g + x T e y J, 

X2=\[x T e- yJ ^ +x T e-y'}. (8) 



In the above equation xt = 2pr/ ' \fs and xt = \J x T + At with r = M^/s. Note that since 
in this process, the J/ip and the photon are both detected, the values of pr, y 3 ^ and 
y 1 are all known and therefore the kinematics is determined completely. The measured 
cross-section is, therefore, fully differential in the kinematic variables. 

We will also discuss the partially integrated cross-section da/dpx, which is obtained 
by integrating the differential cross-section given in eq. (|7|) over restricted regions of phase 
space. For a given value of xt the allowed range of phase space is given by 



\y 



J ^\ < cosh" 1 



x T 



_ h( az^)<^< ta( iz^!L_ ) . (9) 

In order to illustrate the effectiveness of this process in differentiating between different 
models of the EMC effect, we have chosen three representative cases. These are a) the gas 



model JTT|, b) the rescaling model and c) the six-quark cluster model JT^]. As was 
already mentioned, in all these models the construction of nuclear parton densities starts 
from a choice of nucleon parton densities. The model parameters are then fixed, using the 
DIS data with nuclear targets. We refer the reader to the original papers for the details 
of the models and their fits to the DIS data. The model (b) has no cumulative effects 
whereas the others are two-component models where the nuclear structure function has 
two contributions - one from the normal nucleons and the other from nucleons that share 
their partons with the other nucleons. In Fig. 1 we plot the ratio p g defined in eq.(^) 
as a function of X2 for the three models for A = 56[|. It may be noted here, however, 
that p g does not depend very strongly on the mass number A of the nuclear target. All 
the models give somewhat similar predictions, in the vicinity of 1.0, for p g in the range 
0.1 < X2 < 0.2. The deviation of the ratio from unity is significant for all the models at 
larger values of x<± and also at the smaller values for the rescaling model. It should be 
noted here that the choice of nucleon parton densities in each model is different. As a 
result, the relative behaviour of g A in different models can be different from that for p g 
shown here. The axis on the top of the figure shows, as an illustration, values for the 
symmetric configuration y 1 = 0, y J ^ = for the FNAL beam energy of 800 GeV in the 
fixed-target mode, i.e., y/s = 38.75 GeV. As is clear from eq. (0) the ratio of the triple 

*A similar figure appears in the paper by Godbole and Sridhar (Fig. 1 of that paper) quoted in Ref. jjj. 
The curves for the gas model and the rescaling model shown there are incorrect and should be replaced 
with those shown here. 
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differential cross-sections for pA and pp, is directly proportional to p g . Hence this figure 
also shows the expected ratio of the cross-sections as a function of p? (for the fixed-target 
case) for the symmetric point that we have chosen. Moving away from the symmetric 
configuration in the rapidity space to negative values, it is possible to map large x 2 values 
to lower pt values than shown in Fig. 1, where the cross-section is still appreciable. 
Conversely, it is possible to probe the smaller values of x 2 by going to positive rapidities. 
In this way one can obtain information on gluon densities over a fairly large range of x 
values. 

In Fig. 2 we first show the cross-section da/dpx, for the process p + A — > J/^ + 7 — >■ 
Z + /~7 , obtained by integrating eq. (|7|) over the region \y J ^ ,1 \ < y ma x, for the three models 
for the fixed-target mode of FNAL, i.e. yfs = 38.75 GeV, and also for AA collisions for 
RHIC for v/i = 100 GeV. We have chosen y max = 2.0 for FNAL and y max = 4.0 for RHIC. 
Including both electrons and muons will enhance the cross-section by a factor of 2. We 
see that, since the ^-dependence of g A is different in different models, the pr-dependence 
of the cross-sections in these models is also quite different. Recall that the assumed 
nucleon gluon density in different models is different. Therefore the gas model predicts a 
harder px distribution than the quark-cluster model, contrary to what one would expect 
from Fig. 1. Interestingly, even the lowest prediction for the FNAL energy, (that of the 
rescaling model) corresponds to about 30-40 pb at pt = 3 GeV. With the total integrated 
luminosity of 100 pb -1 that is foreseen in this case, this cross-section corresponds to large 
enough statistics to facilitate a good study of the process. Even for pp collisions we should 
have ~ 100 events at pr = 3 GeV. 

The increased value of y/s at RHIC causes the cross-section to be dominated by the 
low x region. The px distribution is much broader, of course. The much higher cross- 
sections in this case compensate for the decrease in luminosity by about a factor of 100 
in going from pp to AA collisions. Since the g A distribution is much broader for the gas 
model, the increase in going from FNAL energies to RHIC is only a modest factor of 5 
(after taking into account the additional factor of A for the RHIC case). For the other 
two models the cross-section increases by more than an order of magnitude. As a matter 
of fact the relative size of the cross-sections predicted in the models changes completely 
as one goes from FNAL to RHIC. At RHIC largest (smallest) cross-sections are foreseen 
for the rescaling (gas) model, which is exactly opposite to the situation at FNAL. Thus 
combining the information on J/ip + 7 production in the fixed target experiments and 
at RHIC, one can probe g A over a wide range of x 2 and can indeed discriminate between 
different models. Of course, owing to the completely differential nature of the measured 
cross-section information can be obtained on p g (x 2 ) over a wide range of x 2 even at a 
fixed value of y/s, by studying the triple differential cross-sections of eq. (|7|). We choose 
Pt = 3 GeV for doing this: at these large values of pr, the prompt photon detection is 
free of background and the rates are still appreciable. The decay leptons in this case will 
also have on the average pr — 2 GeV. 

In Fig. 3 we show the triple differential cross-section for the process for A = 56, for 
the gas model as a function of y J ^ and y 7 , for px = 3 GeV at FNAL energy. We have 
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restricted ourselves to values of y J ^ and y" 1 in experimentally feasible ranges, \y J ^, y 7 | < 
2.0, even though the kinematic limit is somewhat higher. As is clear from the figure, most 
of the cross-section is contained in this region. Indeed, the figure also shows that fairly 
large values of cross-sections are possible at large negative and positive rapidities. These 
curves show a certain asymmetry between y J ^ and y 1 . This is a result of the fact that the 
nuclear gluon density is much harder than the nucleon gluon density in this model. For 
the quark-cluster model, a very similar configuration results. Only the absolute values of 
the cross-section are somewhat lower. The rescaling model predicts the least cross-section 
of the three models. This is a reflection of the softer nuclear gluon density in the rescaling 
model and the somewhat large value of pt chosen here. The softer gluon also results in 
narrower rapidity distributions in this case. Thus one sees that the differences in the g A 
are indeed reflected in these distributions. 

To get a clearer picture of the differences in the predictions between different models, 
as well as to gauge the measurability of the process, we show in Figs. 4(a), (b) and (c) 
the contours of constant triple differential cross-sections (in pb/GeV) in the 

yJ/iP _ yl 

plane, again for = 3 GeV and for the fixed-target case, for the gas model, the rescaling 
model and the quark-cluster model respectively. The asymmetry between y J ^ and y 1 is 
significantly different in all three cases. The gas model predicts measurable cross-sections 
(~ 1 pb) all the way up to \y J ^,y" / \ ~ 1.6, whereas in the rescaling model these are 
obtained only up to rather small rapidities of 0.8 or so. As was said before, this reflects 
the much softer g A in this case. 

Superimposed on these contours are the contours of constant %2- Since different pairs 
of y J ^, y^ correspond to the same 22 at a given pt, the ratios of triple differential cross- 
sections integrated along this contour will yield p g at a fixed value of 22- Note that 
while taking the ratios, factors of X\ and px in eq. (|7]) cancel; we can thus integrate the 
ratio along the Xi contour. This will clearly help increase the statistics and reduce the 
errors without losing the differential nature of the information. This enhancement of the 
statistics is particularly relevant for proton targets, since it is necessary to take the ratios 
of triple differential cross-sections for the pA and pp case to obtain p g . From Fig. 4(a), 
for instance, we see that with the assumed luminosity of 100 pb -1 at FNAL, for the gas 
model we expect ~ 600 events corresponding to X2 = 0.3. With the proton (deuterium) 
target we expect ~ 15(30) events. A further such integration is possible along the pt axis 
as well, although the cross-section falls rather rapidly in pt, as shown in Fig. 2. Note 
that for an effective study of p g it is necessary that we choose the lighter target to be 
proton or deuterium, as the EMC effect saturates very fast with increasing A. 

In conclusion we see that the associated production of J / ip and photon at large px is 
a very sensitive probe of the gluon densities. In pA fixed-target experiments, this process 
can be used to study the gluonic analogue of the EMC effect. As we have shown in this 
letter, a study of the ratio of the triple differential cross-sections, integrated along a fixed 
£2 contour, for pA and pp collisions will directly yield the ratio of gluon distributions 
in the nucleus and the nucleon over a wide range of x. This information by itself can 
be used to discriminate between the different models of the EMC effect. The power of 
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discrimination can be increased further by combining this with the measurement of this 
process at RHIC, which is sensitive to g A at smaller values of x ~ 0.1. In this case it will 
not even be necessary to measure the triple differential cross-section and a measurement 
of da/dpx will suffice to distinguish between the different models. 
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Figure captions 



Fig. 1 The predictions of the gas rescaling [12|, and quark-cluster model f!3 |, for the 
ratio of the nuclear gluon density g A to that of the nucleon p g , for A = 56, as 
a function of x, the momentum fraction. The corresponding values of pr for the 
FNAL fixed-target energy, yfs = 38.75 GeV, are shown on the top axis. 



Fig. 2 The expected transverse momentum spectrum of the J/ip and 7 produced in the 
pA (AA) collisions for the FNAL (RHIC), at y/s = 38.75 (100) GeV, for A = 56, for 
the gas, rescaling and quark-cluster model. The spectrum is obtained by integrating 
over \y J ^,y J \ < 2.0 (4.0) for FNAL(RHIC) and the leptonic branching ratio of the 
J/ij) into one lepton species has been folded in. 



Fig. 3 The triple differential cross-section da/dpr dy^dy 3 ^ (in pb/GeV), expected in the 
gas model, for p + A — > + 7 — ► / + /~7, at the fixed-target FNAL energy 
(y/s = 38.75 GeV), as a function of the rapidities y J ^ and y 1 of the J/ip and 
the 7 respectively, at pt = 3 GeV and A = 56. 



Fig. 4 Contours of constant triple differential cross-section, da/dpr dy^dy -7 ^, for p + A — ► 
J/i/j + 7 -> at v/i = 38.75 GeV, for A = 5Q,p T = 3 GeV, for the gas 

(a), rescaling (b) and quark-cluster model (c) of the nuclear parton densities. The 
values shown are in pb/GeV. Also superimposed are the contours of constant x%, 
the momentum fraction carried by the nuclear partons. 
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